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ABSTRACT 

Solar neutrinos coming from different nuclear reactions are now detected with a high statistics. Con- 
sequently, an accurate spectroscopic analysis of the neutrino fluxes arriving on the Earth's detectors 
become available, in the context of neutrino oscillations. In this work, we explore the possibility of 
using this information to infer the radial profile of the electronic density in the solar core. So, we 
discuss the constraints on the Sun's density and chemical composition that can be determined from 
solar neutrino observations. This approach constitutes an independent and alternative diagnostic to 
the helioseismic investigations already done. The direct inversion method, that we propose to get the 
radial solar electronic density profile, is almost independent of the solar model. 

Subject headings: Neutrinos - Sumevolution -Sundnterior - Stars: evolution -Stars:interiors 
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1. INTRODUCTION 

Neutrinos, once produced in the core of the Sun, reach 
the Earth in less than 8 minutes. On their interplanetary 
journey, solar neutrinos very rarely interact with other 
particles. Luckily for physicists, neutrinos with a scat- 
tering cros s-section varyi ng from some 10~*^cm^ to some 
10~^^cm^ ()Bahcalllll989[ ). occasionally interact with the 
baryons present in the solar neutrino detectors. There- 
fore, the neutrinos produced in the nuclear reactions of 
proton-proton (PP) chains and Carbon-Nitrogen-Oxygen 
(CNO) cycle are natural probe-particles of the constitu- 
tive plasma of the Sun's core. 

The current theory of neutrino physics, following in 
the footsteps of the fundamental theory of elementary 
particles, has shown that neutrinos appear in three (at 
least) flavours (jHernandezli20iCl( ) : electron- neutrino {v^), 
muon-neutrino [vp) and tau-neutrino {vt)- The theory 
states that by vacuum flavour oscillations, when moving 
in vacuum, neutrinos have the ability to switch cycli- 
cally between different flavours. This mechanism of neu- 
trino oscillations was proposed by Pontecorvo around 
1940. In the late seventies, Lincoln Wolfenstein and col- 
leagues suggested another mechanism for neutrinos to 
change flavou r, as neutrinos mo ve through a high den- 
sity medium (jWolfensteinl 119781) : in an environment of 
high density material, the effective mass of the propagat- 
ing neutrinos is different from the mass of neutrinos that 
propagates in a vacuum. Since the oscillations between 
different flavours of neutrinos depend on their mass, 
the oscillations of neutrinos in dense media are differ- 
ent from neutrino oscillations in vacuum. This process is 
now known as Mikheyev-Smirnov- Wolfenstein (MSW) or 
"matter oscillations". These two mechanisms of flavour 
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oscillations are central to the modern theory of neutri- 
nos. Through these mechanisms the electron- neutrinos 
produced in the Sun's core have a non-zero probability of 
being detected as a muon-neutrino or tau-neutrino when 
they reach the Earth. In particular a small fraction of 
electron-neutrinos changes flavour in the Sun's interior 
when they propagate through the high density plasma 
medium of the Sun's core. 

In recent years, significant progress has been achieved 
in the understanding of the m echanism responsi- 
ble for the neutri n o o scillations fEgu chi et al.l 120031 : 
iFukuda et all 120011 : [Ahmad ct al. 2001). The current 
theory of neutrino physics successfully explains the neu- 
trino observations. Furthermore, the theory has been 
used successfully to determine the fu ndamental proper- 
ties of neutrinos (jSchwetz et al.|[2008() . 

If, during several decades, the solar neutrino fluxes ar- 
riving on Earth have been difficult to interpret, this is 
no longer the case as the flavour oscillations at differ- 
ent energies are well understood. In particular, the pro- 
duction of these neutrinos has been confirmed by helio- 
seismic data. One may consider as a success that these 
two disciplines (solar neutri nos and he lioscismology) 
agree within the error bars (see Turck-Chie zc et al. 201(| 
iTurck-Chieze fc Couvida^l201lL and references therein). 



Such an agreement, will now allows us to progress a step 
further. We believe that solar neutrinos will very likely 
become a powerful tool to probe the core of the Sun be- 
cause the production of neutrinos is very sensitive to the 
local properties of the solar plasma. 

In this paper, we focus on the sensitivity of solar neu- 
trinos to the plasma in the solar core, and explore the 
possibilities of probing the physical properties of such 
plasma by using the solar neutrino flux measurements. 
This point has been proposed by John Bahcall and Ray- 
mond Davis about 50 years ago (Davis .1964). In prin- 
ciple, the structure of the solar core can be studied 
by means of neutrino spectroscopy in two fundamen- 
tal ways: (i) to diagnose the temperature profile in the 
Sun's core, by measuring the total number of electron- 
neutrinos produced in each nuclear reaction of PP chains 
and CNO cycle, in particular the boron flux that is the 
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Fig. 1. — The electron-neutrino fluxes produced in the vari- 
ous nuclear reactions of the PP chains and CNO cycle. The 
center of the different neutrino sources ^j(r), {r)j (with j = 
pep,pp,^ BJ Be,^^ N,^^ O,^^ F) are the following ones (in unity of 
Rq): 0.099 (pep), 0.117 (pp),0.125 (Hep), 0.054 (^B), 0.070 CBe), 
0.074 ("AT), 0.055(1^50), 0.053(l''F). These neutrino fluxes were 
calculated for the standard solar model, using the most updated 
microscopic physics data. This solar model is in agreement with the 
most current helioseismology diagnostic and other solar standard 
models published in the literature (see text). For each neutrino 
type j, ^j{r) = (l/Fj) dfj{r)/dr is drawn as a function of the 
fractional radius r for which fj is the flux in and Fj is the 
total flux for this neutrino type. 

most sensitive to the temperature (|Turck-Chieze et alJ 
I2012t iTurck-Chieze fc LopesI I2012D and (ii) to measure 
the radial electronic density of matter of the Sun by de- 
termining the amount of electron-neutrinos that are con- 
verted into another flavour. This article is a first attempt 
to study this second physical information through neu- 
trino probes. 

The neutrino flavour oscillations by the MSW mech- 
anism are particularly significant in the solar interior, 
where there is a wide radial varia tion of the plasma den- 
sity. The standard solar model (|Turck-Chieze fc LopesI 
Il993| ) partly validated by helioseismology, predicts that 
the density inside the Sun varies from about 150 g cm~'^ 
in the centre of the star, to 1 g cm~^ at half of the solar 
radius. 

What makes this diagnostic particularly powerful is 
the possibility of getting a direct measurement of the 
plasma electronic density almost independently of the 
solar model. The strong dependence of neutrino oscil- 
lations with the local electronic density of matter opens 
the possibility of inferring the matter density and chemi- 
cal composition in the nuclear region. Since the neutrino 
oscillations in matter depend strongly on the matter den- 
sity and weakly on the chemical composition, this is at 
first, a diagnosis of the radial density profile of matter 
in the solar core. The composition diagnostic can come 
from precise CNO neutrino detections. Therefore, we can 
anticipate that with the expected increased accuracy of 
the solar neutrino detections in the future, it will be pos- 
sible to extend this analysis to the solar core chemical 
composition determination. 

2. CURRENT MODEL OF NEUTRINO PHYSICS 
OSCILLATIONS 



T he theory of neutr i no oscillations in vac uum and mat- 
ter ()Hernandezl l2010t iSchwetz et al.l I2008D has success- 
fully addressed the problem of the solar neutrino deficit 
- the discrepancy between the neutrino flux detected 
and the theoretic al predictions of t he solar standard 
model (iTurck-Chieze fc LopesI [Tool iFogh et al.l [200l 
iGonzalez-Garcia k Maltonil |2008[) . The theory pro- 
vides a theoretical solution in full agreement w i th all 
solar neutrino exper iments ()Bellini et alJ l2011bl . I201CII : 
lAharmim et al.ir201Clt ) as well as with the data obtained 
by t he KamLAND reactor experiment (Eguchi ct ^ 
120031 ). Presently, the theoretical model for neutrino 
flavour oscillations is defined by means of six parame- 
ters: the difference of the squared characteristic masses 
{Aml2, ^mli), the mixing angles (^12, ^23, ^31) and 
the CP-violation phase. The mass square differences 
and mixing angle s are known with reasonable accuracy 
(jFogli et al.ll2009| ): Am|]^ is obtained from the experi- 
ments of atmospheric neutrinos and Amjj is obtained 
from solar neutrino experiments. 

The mixing angles are not uniformly well defined: 812 is 
obtained from solar neutrino experiments with an excel- 
lent precision; ^23 is obtained from the atmospheric neu- 
trino experiments, this is the mixing angle of the highest 
value; has been first estimat ed from the Chooz re- 
actor (jApollonio fc Baldinilll999D. its value is v ery small 
and was still very unccrtain jjFogh et al.ll2009D . Nowa- 
days with Daya Bay and Reno, the situation is largely 
improved (Fogli et al. 2012), as shown in the next para- 
graph. However, present expe riments cannot fi x the 
value of the CP- violation phase (jHernandedl2010( ). 

An overall fit to the data obtained from the dif- 
ferent neutrino experiments: solar neutrino detectors, 
accelerators, atmospheric neutrino detectors and nu- 
clear reactor experiments suggests that the param- 
eters of neutrino oscillations are the following ones 
(|Gonzalez-Garcia fc Maltonil I2008D : Ato|i ~ 2.46 ± 
0.12 Amfa ~ 7.59 ± 0.2010-5ey2, 612 = 

(34.4 ± 1)°, 023 = (42.8t^;^)° and ^13 = (5.6t^-0)°. The 
recent progress leads to ^13 = (8.6^Q;4g)° and Amgj ~ 

-2.43+g:;^^10^^ey^ (|Gonzalez-GarciaeraII[2?)Tl . 

In the limiting case where the value of the mass dif- 
ferences, Ato^2 or Ato§]^ is large, or one of the angles 
of mixing (^12, 6*23, 6*31) is small, the theory of three neu- 
trino fiavour oscillations reverts to an effective theory of 
two neutrino fiavour oscillations (Hernandez 2010). Bal- 
antekin and Yuksel have shown that the survival proba- 
bility of solar neutrinos calculated in a model with two 
neutrino flavour oscillations or thre e neutrino flavour os- 
cillat ions have very close values (jBalantekin fc Yuksel! 

Hool. 

In the present work, for reasons of convenience and 
simplicity, we will restrict our study to the theory of two 
neutrino flavour oscillations. The generalization of the 
results to a theory of three neutrino flavour oscillation s 
is obvious ()Lisi fc Montaninol 119971: iGando et all 120 111) . 
The survival probability of electron-neutrino function 
P^^(= P{iye — > t/e) ) in a two flavour neutrino theory 
is given by 

Pu, {E, r) = ^+ F(7) cos (20„) cos (20„O (1) 
where Oy is the angle of oscillation in vacuum and 9m 
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is the angle oscillation in matter. The function ^(7) 
is a first-order correction to the adiabatic approxima- 
tion of the propagation of neutrinos in matter, where 
7 is the adiabatic parameter which depends of the lo- 
cal properties of the plasma (jLandau fc RosenkewitschI 
[1932) ■ The adiabaticity occurs when the electronic den- 
sity of the propagating medium is a slow varying func- 
tion over the neutrino travel path (or equivalently the 
solar radius). This approximation is valid in the solar 
interior, in the case where 7 S> 1, a condition that is 
verifi ed in much of the solar core and the radmtive re- 
gion (jLandau fc RosenkewitschI 119321 : iZeneri 1193^ . As 
F(7) = 0.5 — P12 where P12 corresponds to the transi- 
tion between two distinct neutrino flavours, at first order 
it is reasonable to consider that P12 = 0. The numeri- 
cal va lue of -^(7) is approximately equal to 0.5 (jBilenkvl 
[2010). The angle 6y depends on the distance from the 
Sun to the Earth. In this model we assume 9^ to be equal 
to 012 ■ The phase 9m is the most important term in this 
analysis, since it depends on the properties of the plasma 
in the Sun's core, namely, the local electron density. The 
neutrino mixing angle 9m is given by 



sin 29,n 



sin 29„ 



^ {V/Am^ - cos (26'„))2 + sin^ (29^) 



(2) 



where Am^ is the difference of the squared masses we 
consider to be equal to Am^2- The function of the solar 
plasma V is given by 



V{E,r) = 2^/2G'/ ne{r)E 



(3) 



where G/ is the Fermi constant, rie(r) is the electron 
density of plasma and E the energy of the neutrino. The 
electron density ne{r) = No p{r)/^e{r) where /Ze in the 
mean molecular weight per electron, p{r) the density of 
matter and No the Avogadro's number. 

3. NEUTRINO PRODUCTION IN THE SUN'S CORE 

The neutrino fluxes produced in the various nu- 
clear reactions of the PP chains and CNO cycle have 
been computed for an updated version of the solar 
standard model (iTurck-Chieze et al.ll2004t iBahcall et al.l 



2005a\'Guzik fc Mussack'l2 010l:lTurck-Chieze et al.l l2010l: 
Turck-Chiczc fc Couvidat 2 011[ ). Figure [1] shows the lo- 
cation of the different neutrino emission regions of the 
nuclear reactions. In the Sun's core, the neutrino emis- 
sion regions occur in a sequence of shells, following closely 
the location of nuclear reactions, orderly arranged in a 
sequence dependent on their temperature. 

The first reaction of the PP chains, the pp nuclear reac- 
tion, has the largest neutrino emission shell. This region 
extends from the center of the Sun up to 30 % of the 
solar radius. The pep reaction has a neutrino emission 
shell identical to the pp reaction, although only up to 25 
% of the solar radius. These two key nuclear reactions 
are strongly dependent on the total luminosity of the 
star. This is the reason why different solar models with 
the same total luminosity produced the same pp and pep 
neutrino emission shells. The neutrino emission shells of 
^B-v and Be-v extend up to 15 % and 22 % of the solar 
radius. Finally, it is worth noticing that the maximum 
emission of neutrinos for the PP chains nuclear reactions, 
follows an ordered sequence (see figure[T]): ^B-v, Be-v, 
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Fig. 2. — The solar neutrino energy spectrum predicted by the 
standard solar model. The solid curves correspond to the to- 
tal neutrino fluxes produced in the various nuclear reactions of 
the PP chain and CNO cycle. The dashed curves correspond to 
electron-neutrino fluxes of the various nuclear reactions after neu- 
trino flavour conversion. The neutrino fluxes from continuum nu- 
clear sources are given in units of cm~^s~^Mev~^ . The line fluxes 
(indicated in the legend with (*)) are given in cm^-^s"-"^. 

pep-v and pp-v with the maximum emission located at 5 
%, 6 %, 8 % and 10 % of the solar radius, respectively. 

The CNO cycle nuclear reactions produce the following 
electron-neutrinos sub-species: ^^O-v, ^'^F-i/ and ^^N-v 
released in emission shells identical to the ^B-v. The 
^^N-i/ have a second emission shell located between 12% 
and 25 % of the Sun's radius. The emission of neutrinos 
for ^^O-v, ^'^F-iy and ^^N-v is maximal at 4%-5 % of the 
solar radius. The ^^N-i/ neutrinos have a second emission 
maximum which is located at 16 % of the solar radius. 

4. NEUTRINO FLAVOUR OSCILLATION IN THE SUN 

The neutrino emission reactions of the PP chains and 
the CNO cycle are produced at high temperatures in dis- 
tinct layers in the Sun's core. Similarly, the neutrino 
flavour oscillations occur in the same regions. The av- 
erage survival probability of electron-neutrinos in each 
nuclear reaction region is given by 



a: 



P^^{E,r)^j{r)ATrp{r)r^dr (4) 



where Aj is a normalization constant given by Aj = 
J ^j{r)ATTp{r)r^ dr and $j (r) is the electron-neutrino 
emission function for the j nuclear reaction, j corre- 
sponds to the following electron-neutrino reaction sub- 
species: pp, pep, ^B, ''Be, ^^N, ^^O and ^''F. 

Figure [2] shows the energetic neutrino spectra (jBahcalll 
[l989) for an updated version of our standard solar model 
before and after the flavour conversion. Figure [3^ shows 
the survival probability of electron-neutrinos produced in 
the regions where occur the different nuclear reactions. 
This survival probability of electron-neutrinos (P^^ {E))j 
is very similar for low- and high-energy neutrinos but 
presents some differences between. Unfortunately, as is 
well known, the emitted neutrino energy spectrum is lim- 
ited to a specific energy range for each nuclear reaction. 
Nevertheless, to highlight the sensitivity of neutrino to 
MSW flavour oscillations, we choose to represent the sur- 
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Fig. 3. — The survival probability of electron- neutrinos in func- 
tion of the neutrino energy for the different regions of emission, a) 
survival probability of electron-neutrinos produced in the different 
regions b) difference between the survival probability of the differ- 
ent electron-neutrinos and the reference curve corresponding the 
central emission. The coloured part of the curves (a), and sim- 
ilarly the bold part of the curves (b), indicate the energy range 
of neutrinos produced in the Sun's core for each nuclear reaction 
(Cf. Figure |2]|. This reference curve defines the survival proba- 
bility of electron-neutrinos in the centre of the Sun for which the 
MSW flavour oscillation mechanism is maximum (dashed red curve 
of figure 3a) . 

vival probability of electron- neutrinos in all the available 
energy range, so that the regions of interest are indicated 
in colour, or by a single colour square in the case of line 
fluxes (Cf. Figure 2). 

The neutrino change of flavour for low-energy neutri- 
nos is due to vacuum oscillations, and for high-energy 
neutrinos it is caused by a cumulative effect of oscillations 
in vacuum and matter (MSW effect). The experimen- 
tal neutrino flux measurements for low-ene rgy neutrinos 
have been used to determine the value of 0i2 ([Bellini et alj 
120101) . The survival probability of electron- neutrinos 
{P^^{E))j with intermediate values of energy, between 
1 MeV and 10 MeV, have a strong dependence on the 
electron density of the plasma. Since the production of 
electron-neutrinos occurs in various nuclear reactions at 
different layers in the Sun's interior, a sharp differentia- 
tion is observed between the different {P,y^ {E))j curves 
for the neutrinos with intermediate values of energy. Fig- 
ureint shows a well-ordered sequence of curves that corre- 
sponds to the difference between the survival probability 
of electron-neutrinos produced by each nuclear reaction j 



and the survival probability of electron-neutrinos in the 
center of the Sun. This sequence of curves occurs as a 
result of the regular decrease of matter density from the 
center of the Sun. As a consequence the difference be- 
tween the survival probability of electron-neutrinos pro- 
duced by each nuclear reaction j and the survival proba- 
bility of electron-neutrinos produced in the center of the 
star increases with the distance of the location of the 
nuclear reactions to the center (Cf. Figure 3b). For ex- 
ample, in the case of pp electron-neutrinos this difference 
is of the order of 0.06 which is three times larger than in 
the case of electron-neutrinos. This is due to the fact 
that the electron-neutrino source is located near the 
center of the star (Cf. Figure 1). Nevertheless, the pos- 
sibility to observe such effect is somewhat limited, once 
the neutrino flavour oscillations induced by matter de- 
pends on the location of the neutrino source in the Sun's 
core (cf. Figure 1), as well as the energy of the emitted 
electron- neutrinos, i.e., the neutrino emission spectrum 
(cf. Figure 2). In the case of the pp nuclear reaction, 
the emitted neutrinos have a maximum energy of 0.42 
MeV, consequently the survival probability difference is 
smaller than 0.005. This effect is more expressive in the 
case of the ^B nuclear reaction for which the emitted 
neutrinos have a maximum energy of 14.02 MeV. In this 
case, the effect is maximum for neutrinos with a energy 
of 4 MeV which have a survival probability difference of 
0.02. Identical behaviour will be observed for the CNO 
nuclear reactions (namely the nuclear reactions related 
with the production of ^^N , ^^O and chemical ele- 
ments) which emitted neutrinos with maximum energies 
slightly above 1 MeV, for which the survival probability 
difference is of the order of 0.014. An equally pronounced 
difference is observed for the spectral line of pep neutri- 
nos. In the case of the Be neutrino the effect is impor- 
tant for the high-energy neutrino line and very small in 
the case of the small-energy neutrino line. In particular, 
it should be possible to check experimentally the mod- 
ification of the ^B energy neutrino spectrum caused by 
matter fiavour oscillations due to the radial dependence 
of ^B electron-neutrino source. 

5. THE PHYSICS OF THE SUN'S CORE AS PROBED BY 
NEUTRINOS 

The solar model has been checked in most of 
the solar radiative region by means of the high 
pre cision data of SOHO he l ioseismic instruments 
(Tu rck-Chieze fc Couvida^ [20TTI : iTurck-Chieze fc LopesI 
I2012i) . 

As shown in equations 1 to 3, the radial profile of the 
electron density of the solar standard model is a funda- 
mental ingredient to test the neutrino physics theory. Up 
to now this quantity has been checked by the detection 
of the first gravity modes wh ich are really sensitive to 
the central region of t he Su n ()Turck-Chieze et al.ll20T2l: 
iTurck-Chieze k, LoDesl[20T^ . In the next few years, with 
the increase of accuracy on the measurements of several 
solar neutrino experiments, the situation could reverse: 
neutrino fiuxes should start to be used to diagnose and 
to infer the thermodynamic properties of the Sun's core. 

The neutrino fiux emissions inside the Sun are sensitive 
to the local values of temperature (specifically some of 
them), matter density, chemical composition and elec- 
tronic density where the nuclear reactions are taking 
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Central parameters of standard and modified solar models 



Solar Models'' SSM Sun A Sun B 



Model Values 

density (g/cm^) 151 161 167 

mean molecular 

weight per electron fj,e 1.69 1.65 1.64 



Neutrino Energy (MeV) 

Lower Cut-off Ert> 2.0 1.83 1.75 



Inverted Values 

density (g/cm'^) 

with /ie solar model 151.5 [0.3%] 161 [+0.4%] 168 [+0.6%] 
density (g/cm'') 

with /ie = /isSM = 1-69 151.5 [0.3%] 165 [+2.8%] 173 [+3.6%] 



a 

The solar standard model (SSM) and the modified models, Sun A 
and Sun B, were calibrated in order to keep the total luminosity of 
the present Sun ( Figure|4]shows the corresponding density profile). 

^ Ere is the minimum energy that a neutrino must have to be 
affected by the MSW flavour oscillation mechanism. The values 
between [. . .] show the percentage difference between the model 
values and the inverted values. 

place. 

In particular, as the IMSW flavour oscillation is 
strongly dependent on the local electronic density, it is 
possible to infer this quantity from the survival probabil- 
ity of electron- neutrinos. Furthermore, since electronic 
density depends on matter density and chemical compo- 
sition, it will also be possible to obtain some information 
about these quantities. 

5.1. The sensitivity of neutrino flavour oscillations to 
the central density 

Any standard solar model calibrated for the present 
day total solar luminosity has a radial profile of temper- 
ature and density in the core that results from the bal- 
ance found between the energy produced by the nuclear 
reactions and the energy transported to the surface. So- 
lar models with slightly different physical assumptions, 
have different radial profiles of temperature and matter 
density among other quantities. This is the result of the 
readjustment of the Sun's internal structure caused by 
the need to obtain the same total luminosity. 

Our evolution code is an up-to-date version of the 
one-d imensional stellar evolution code CESAJM (jlVlorell 
I1997D . The code has an up-to-date and very refined 
microscopic physics (updated equation of state, opaci- 
ties, nuclear reactions rates, and an accurate treatment 
of microscopic diffus ion of heavy e l ements), in cluding 
the solar mixture of lAsplund et all (|200a . I2009[ ). The 
solar models are calibrated to the present solar radius 
Rq = 6.9599 X 10^° cm, luminosity Lq = 3.846 x 
Iq33 ergs~^, mass M© = 1.989 x 10 ^^ g, and age 
to, = 4.54 ± 0.04 Gvr (e.g. [l\ir ck-Ch ieze fc CouvidatI 
[20TltlTurck-Chieze. Piau fc Couvidatll2nilH) . The mod- 
els are required to have a fixed value of the photo- 
spheric ratio {Z/X)q, where X and Z are the mass 
fraction of hydrogen and the mass fraction of elements 
heavier than helium. The value of {Z/X)q is de- 
termined acco r ding to the solar mixture proposed by 
lAsplund et al.l (|2005[) . Our reference model is a solar 



standard model that shows acoustic seismic diagnostics 
and sol ar neutrino fluxes near from other solar standard 
models fTurck-Chieze & Lopes"1993VTurck-Chieze et aU 
2004; Bahcall et al. 2005b; Guzik & IViussack 20X3 
iSerenelh et al.ll2009t iTurck-Chic ze et al.1 12010'). 

In Table 1 we present the central density and the mean 
electronic molecular weight of this standard solar model 
and of two other non-standard solar models. These last 
two models have their radiative energy transport slightly 
modified in the core in order to obtain a Sun's model 
with a higher central density. Figure |4^ shows the ra- 
dial density profile in the Sun's core for the three mod- 
els. The amount of neutrinos converted to the non- 
electron flavour by the IVISW flavour oscillation mech- 
anism is dependent on the central density profile. Fig- 
ure IIId shows the difference between the survival prob- 
ability of electron-neutrinos for each of the two modi- 
fied solar models relatively to the standard solar model. 
The difference is more significant for the model with the 
highest central density. It is clearly illustrated that an 
increase of the central density leads to an increase in the 
neutrinos converted by the IVISW flavour oscillation in 
all electron-neutrino subspecies. 

The survival probability of electron-neutrinos can be 
used to infer the matter density in the centre of the Sun, 
once that precise neutrino measurements will become 
available. In the following we describe a simple proce- 
dure that highlights the sensitivity of neutrino flavour 
oscillations to the Sun's core density. Only neutrinos 
with energy above a given minimum value will have their 
flavour changed by the IMSW flavour oscillation mecha- 
nism. Several resonances occur in equation ^ for neutri- 
nos with energy Er such that V{Er, r) — Arri^ cos (20^,). 
For the lowest neutrino energy value that verifles the 
previous equation, we choose to call it E,-c- the cut-off 
minimum resonance neutrino energy for a given solar 
model. The value of E,- is minimum for the maximum 
value of the density of matter (equation [3]) . Therefore, 
each solar model has a unique cut-off minimum neutrino 
energy E^c given by 

_ Ato2cos(26I^) 

-t^TC ^ ( O ) 

2V2GfNoPc/fiec 

where pc in the central matter density and fiec is the 
central value of mean molecular weight per electron. 

In Table 1, we show the value of E^^ for the three solar 
models computed from the survival probability electron- 
neutrino function: if the resonance condition is veri- 
fied, 9m = 7r/4 and the survival probability of electron- 
neutrinos (equation [ij has the specific value of 0.5. The 
value Ere is determined for each solar model by identi- 
fying the neutrino minimum energy value for which the 
survival probability of electron-neutrino is equal to 0.5. 
The models with the higher central densities have the 
lowest values of E^c (see Table 1). The central density 
can be calculated from equation ([S]) once the value of Ej-^ 
is known for each solar model and assuming that fiec (or 
the chemical composition) is known. In Table 1 we show 
the "inferred values" of the central density. The small 
variation /.tec between different models allows us to make 
a reasonable estimation of the central density, assuming 
that the variation of resonance energy Ere is due solely 
to variations in the central density. The central density is 
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Fig. 4. — The standard solar model and two modified solar models 
(see table 1): (a) density as a function of radius in the solar core; 
(b) difference of the survival probability of electron-neutrinos of 
two modified solar models and the survival probability of electron- 
neutrinos of the standard solar model. The survival probability of 
electron-neutrinos was computed for neutrinos propagating in the 
center of the star. The continuous-black curve corresponds to the 
standard solar model, the dashed-red curve and dotted-blue curve 
correspond to the solar models, Sun A and Sun B (see table 1). 
In the figure, the vertical lines indicate the location of the present 
solar neutrino measurements: ^ Be-f (0.56 MeV), pep-v (1.44 MeV) 
and ^B-u (8.9 MeV). 

inverted at most with a difference of a few percent higher 
than the value of the model, even in the case that is 
poorly known. 

In fact, the central value of the density in the Sun 
might be more difficult to obtain since the electron- 
neutrinos are not produced just in the center, and the 
different ranges of energy of neutrinos are limited (Cf. 
Figure An interesting possibility is to probe the solar 
core by using the neutrinos produced near the centre of 
the star, such as ^B-v, ^^O-v and F-v neutrino fluxes. 
The simultaneous measurement of two different neutrino 
fluxes produced in the same layer of the Sun's interior, 
such as ^B-v and ^^O-u, could be used to estimate simul- 
taneously the matter density and the average molecular 
weight per electron in the core of the Sun. 

5.2. The sensitivity of neutrino flavour oscillations to 
the radial density profile 



All the solar models validated by acoustic and grav- 
ity mode detections, predict that the density of matter 
decreases by 90% from its central value to the value at 
30% R0. The neutrino fluxes produced by the different 
nuclear reactions are sensitive to the local values of the 
electronic density. The different neutrino fluxes '^jir) 
can be interpreted as a local average of nj (r 

In principle taking into account the high sensitivity 
of solar neutrino fluxes to the temperature changes in 
the solar core, we should expect important variations in 
the location of the neutrino source in the Sun's core. 
However, the effect is actually quite small. The rea- 
son is related with the calibration procedure that so- 
lar models are subjected to, which require that the Sun 
in the present age must have the observed solar lumi- 
nosity. Consequently, the production of neutrinos inside 
the Sun, which are strongly dependent on the luminos- 
ity of the star, namely, the neutrinos produced in the 
proton-proton chain nuclear reactions, such as ^B and 
Be, occur sensibly in the same location, i.e., at the 
same distance from the center of the star. As an ex- 
ample, solar models which experiment a decrease of 10% 
in the central density (and a 3% reduction in the cen- 
tral temperature) have a variation of the maximum of 
any of the neutrino sources (both proton-proton chain 
nuclear reactions and CNO nuclear reactions) displaced 
by an amount smaller than 1% of the solar radius, al- 
though the neutrino fluxes change significantly due to 
their high sensitivity to the temperature. Actually, the 
proton-proton neutrino source has the same location, as 
this nuclear reaction, more than any other depends di- 
rectly from the total luminosity of the star. For the pur- 
pose of this article the neutrino sources of different so- 
lar models are considered to be located within the same 
fraction of the solar radius. Therefore, if the funda- 
mental parameters of neutrino oscillations are known, 
the probability of survival of electron neutrinos can be 
used to measure the radial electronic density profile in 
the solar core, i.e., the matter density and molecular 
weight per electron in several layers of the Sun. A sim- 
ple procedure to obtain the electronic density is now 
proposed. We consider first an electronic-neutrino sur- 
vival probability pj determined, at a specific neutrino 
energy Ej, by some neutrino experiment. One assumes 
that pj is determined fully independently from the solar 
standard model by a method identical to the one de- 
scribed in Ber ezinskv fc Lissial (|2b01.) . Preferentially, Ej 
is the value chosen within the energy interval for which 
the MSW neutrino flavour oscillations are significant (cf. 
Figure [3]). For each pair {Ej,pj), the value of the elec- 
tronic density for a certain layer can be computed from 
equations ([T]|3]). We estimate fij = {ne{r))j correspond- 
ing to electronic density of several layers in the Sun's 
core, as fij = l/{2y/2Gf) x aj{V)^/Ej, where {V) ^ 
is computed from pj using the previous equations, and 
aj is a weight-correction parameter unique for each j- 
nuclear reaction, aj is estimated for the solar standard 
model and it takes a value of the order of a unit, in- 
creasing slightly when the neutrino source $j(r) moves 

^ nj(r) is the estimated mean value of the electronic density 
for the layer at the radius rj, where the neutrino source 4"j(r) is 
maximum (Cf. figure [TJ. 



7 



away from the centre of the star. Siniilarly, by using 
the experimental value pj, which is different from the 
theoretical prediction pj, we estimate the value Ti j as a 
correction to the theoretical value rij, from the equa- 
tion Anj/rij — Pj l^pj/pj, where An^ — fij — nj, 
Apj = Pj — Pj and /3j is a coefficient computed from 
the solar standard model. The previous expression is ob- 
tained from a perturbation analysis of equations ([T][3]). 
In this analysis Uj is the inverted value deduced from 
the experimental data. 

Figure [5] shows the inverted electronic density val- 
ues obtained from the survival probability of electron- 
neutrinos as described previously. There is a good agree- 
ment between the values of the electronic density ob- 
tained by inversion and the electronic density of the stan- 
dard solar model. On Figure [5^ we show the electronic 
density values that one can deduce from the survival 
probabilities obtained from the neutrino detector mea- 
surements. More specifically, we present the electronic 
density inverted f rom the Be neutri no fluxes using the 
measurements of iBellini et al.l ()201C1[) a nd the neu- 
trino fiuxes using the measurements of lAharmim et aP 
dlQQT). 

The Borexino experiment measures ^ Be solar neutrino 
rates with an accuracy better than 5%. This corre- 
sponds to a '^Be neutrino flux ^(^Be) = 'i-SltoH x 
10^ cm~^s~^, under the assumption of the MSW- 
LMA s cenario of solar neutrin o oscillations (|Bellini et al.l 
I2011at lArpesella et al.l 120081 ). The estimated survival 
probability, assuming a high metallicity solar standard 
model, was initially estimated to be P^, ( -^^) = 0- 56ib0.1 
at the energy 0.862 McV ( ArpcscUa'ct al.''2008'). Lat- 
ter the Borexino coUaborationf Bellini et al. 2011a) up- 
dated their estimation to P^^^Be) = 0.52 t °;°^. Sim- 
ilarly, as proposed by iBer ezinsky fc Lissial (|2001| ). we 
compute the survival probability for ^B as P,j^{^B) = 
^i^A^B)/^ti^B), where ^t^B) is the total neutrino 
flux integrating all neutrino flavours measured by 
SNO (lAharmim et al.l[2010l) : $t(*B) = 5.046j:J] p^ x 
10^ cm~^s~^ and is the electron-neutrino flri x 

measured by the same detector (jAharmim et al.ll2007[ ): 
$^^(8^) = l-76l[];}^ X 10^ cm-2s-i. It follows that 
P^^ (8B) = 0.35±0.047 for SNO. It is clear that presently 
the error bars are still too large to really estimate the 
small gap with the standard model values. But one 
may hope a progress on the future as we have already 
for ^B neutrinos several improved detections: the neu- 
trino flux measured by the Kamiokande-III experiment 
(|Abe et al.ll20Tl . $^,(^-8) = 2.32j:°:°^ x 10^ cm'^s"!, 
the neutrino flux measured by the Borexino experiment 
is = 2.4t° j X 10^ cm-2s-i (IBellim et al.ll2010n . 

Moreover, we already measure the energy dependence of 
these fluxes down to 3 MeV. So, to illustrate the po- 
tential of this diagnostic and to indicate some research 
perspectives, we present in Figure the electronic den- 
sity values deduced assuming an error bar on the survival 
probability of the order of 4%. It follows that the error 
bar on the electronic density that one can deduce for 
'^Be and ^B fiuxes is largely reduced. This is due to 
the high sensitivity of neutrino flavour oscillations to the 
electronic density. This first study shows the potentiality 
of such type of neutrino diagnostic. Another possibility 
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Fig . 5 . — Radial profile of the electronic density nc (r) (black con- 
tinuous curve) in the core of the Sun of the solar standard model, 
on the same plot appear normalized values of the radial profile of 
the density of matter, central p = 167 g cm~^ (blue dashed curve) 
and the radial profile of molecular weight per electron tJ.e{r), cen- 
tral value of 1.69 (red dotted curve). The points shown in the 
figure correspond to the values of density inverted from the sur- 
vival probability of electron-neutrinos, at the following energies: 
0.56 MeV ( ''Be), 1.44 MeV ( pep), 8.9 MeV (^B), 1.0 MeV C^N), 
1.0 MeV ( i^o) 1.0 MeV {^"^F). The location of the source are 
the following ones: 0.070 Rq C Be, square- green) , 0.099 Rq {PeP, 
square-blue), 0.054 Rq (*_B, square-green), 0.074 Rq {^^N, square- 
magenta), 0.055 -R0 ( ^^O, square-cyan) , 0.053 Rq (^'-F, square- 
red). The inverted electronic density computed from measured ^B 
and ^ Be neutrino fluxes: (a) the values computed from the survival 
neutrino probability obtained from the present results of SNO {^B) 
and Borexino (J Be), (b) idem considering that the experimental 
error in the determination of the survival probability of electron 
neutrinos is 4 %. The red and in cyan circles correspond to the ^B 
and ^ Be neutrino measurements respectively. 

is to fix values of the electron-neutrino survival proba- 
bility curve for neutrinos of low and high energy, as such 
neutrinos have respectively pure vacuum oscillations and 
vacuum plus MSW oscillations. In this case neutrinos 
with energy in the interval 0.5 — 10 MeV, can probe the 
electronic density of different solar layers. 

6. DISCUSSION AND CONCLUSION 

The capability to make such a study is real because 
new neutrino detection experiments are arriving and oth- 
ers are planned for the near future. The diversity of 
neutrino fiux measurements expected at different energy 
levels will allow us to make reliable and detailed neutrino 
inversions. 

In this first study we have used only two mean sur- 
vival probabilities, without taking into account the en- 
ergy dependence of the boron flux and the detailed radial 
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distribution. But the Borexino detector also measures 
the pep neutrino flux, $(pep) = l-6+o 3 ^ 10^cm~^s~^ 
(iBehini et al.l 120121 ) . This preUniinary result is still in- 
sufficiently accurate, but it is not in contradiction with 
the fact that the Sun is hotter than what is suggested 
by the present standard solar model. The pep neutrino 
flux is strongly dependent on the luminosity of the star. 
Therefore, it is an indirect measurement of the total en- 
ergy (luminosity) produced in the nuclear region. So it 
represents a powerful probe to the physics of the nu- 
clear region of the Sun, in para flel to the information 
introduced in the seismic model (jTurck-Chieze fc LopesI 
[Ml . 

Today, there is some uncertainty about the chemical 
composition in the Sun's core. In parallel, all the neu- 
trino measurements agree with the predicted values of 
the sola r seismic model but not so well w ith the standard 
model (jTurck-Chieze &: Couvidatj l2011h . showing some 
capability of understanding the difference with the stan- 
dard model if one extracts the electronic density profile. 
The results obtained so far suggest that in the near fu- 
ture we will be able to probe the physics of the Sun's 
core using neutrinos. 

In this work we have proposed a strategy that allows 
to use the solar neutrino flux measurements to invert the 
electron density of the solar plasma at specific locations 
of the Sun's core. Furthermore, in the near future, it will 
also be possible to put important constraints on the mat- 
ter density and on the molecular weight per electron of 
the solar plasma. The inversion is made based upon the 
assumption that the future Earth neutrino experiments 
(e.g., reactor experiments, superbeams, beta beams and 
neutrino factories) will allow the precise determination 
of basic neutrino oscillation parameters. 

Today, the neutrino fluxes of Be-v, -v and pp-v 
are already well measured and obtained separately. The 
pep-v neutrino flux is also measured although with much 
less precision. It will be interesting to verify if such re- 
sults hold with the increase of accuracy in the observa- 
tions within the next few years. The present electronic- 
neutrino survival probability curve is fixed at low neu- 
trino energy by pp-v, sensitive to pure vacuum oscilla- 
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tions, and by ^B-v, sensitive to vacuum plus MSW os- 
cillations. Therefore it will be relevant to verify if future 
measurements of Be-v and pep-v can give us some infor- 
mation about the electronic density in two distinct layers 
of the Sun, namely 0.1 Rq {pep-v) and 0.07 Rq Be-v). 
There is a real possibility that Borexino or SNO -f exper- 
iments will be able to measure accurately these different 
neutrino fluxes, as well as the neutrino fluxes of some of 
the CNO cycle neutrino emission reactions. 

The LENA (Low Energy Neutrino Astrophysics) solar 
neutrino detector is expected to start to operate within 
the next few years. This detector will be able to perform 
very accurate measurements of the different sources of 
the solar neutrino spectrum, allowing not only very pre- 
cise measurements of the solar core plasma, but also to 
identify possible seasonal neutrino variations. 

In conclusion, we have shown that in the near future 
one may hope to use neutrino spectroscopic measure- 
ments to infer the electronic density of the plasma in the 
core of the Sun. This will be an important and totally in- 
dependent test of the nuclear region of the Sun properties 
in parallel with the helioseismic inversions of matter den- 
sity and sound speed. The proposed method, that could 
be improved, allows to infer the properties of the plasma 
in the very central region of the Sun, which until now 
have only been explored by some gravity mode detec- 
tions with GOLF/SoHO. We have shown in table 1 that 
the inversion of different density profile can be done with 
a reasonable quality if the neutrino fluxes are obtained 
with a good accuracy. Of course this method supposes 
that the neutrino oscillations description is independent 
of solar model predictions. If that becomes the case, one 
will be able to extract from the different solar neutrino 
detections, the electronic density and matter density pro- 
flies with some hints on the CNO composition in the core 
of the Sun. 
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